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Abstract. We give a precise measure of the rate at which repeated differentiation of a ran- 
dom trigonometric polynomial causes the roots of the function to approach equal spacing. 
This can be viewed as a toy model of crystallization in one dimension. In particular we de- 
, termine the asymptotics of the distribution of the roots around the crystalline configuration 

' and find that the distribution is not Gaussian. 

(N 

1. Introduction 

The critical points of an analytic function have a variety of interesting physical inter- 
im-- pretations. For polynomials we have the Gauss electrostatic model: at each zero of the 
polynomial f(z) place identical point charges obeying an inverse linear law. Then the zeros 
of the derivative f'(z) are the points where the field vanishes. To see when this works, just 
write f(z) in factored form and consider the logarithmic derivative f'(z)/f(z). 
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The Gauss model extends to entire functions of order 1, provided one incorporates a 
background field coming from the exponential factors of the Hadamard factorization of the 
function. For such functions there is a general phenomenon that differentiation smooths out 
irregularities in the distribution of zeros. See [S] for details. If the zeros are located in a 
strip around the real axis and their initial distribution is not too irregular, then repeated 
differentiation leads the zeros to approach equal spacing. This regular spacing is known as 
the crystalline configuration [2] and we view the process of repeated differentiation as a toy 
model of crystallization in one dimension. According to the Gauss model, regular spacing 
^ ■ is the equilibrium position and differentiation moves the function toward equilibrium. See 
Figure n for an illustration. 

We now describe the functions we study and then discuss our results. 
A random trigonometric polynomial of degree N is a function of the form 
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N 



'1.1) F(x) = ^^a n cos(nx) + b„ sin(nx) 



n=0 



where a n and b n are random variables. In this paper we will assume that the a n and b n are 
independent real Gaussian distributed with mean and variance a n , and usually we further 
assume that the variances a n are all equal. 

Our concern is with the properties of the real zeros of F(x). Generally F(x) will not 
have all of its zeros real, but the high derivatives of F(x) will have mostly real zeros, and 
those zeros will be close to equally spaced. This can be seen in Figure ^ This is a general 
property |Hj of real entire functions of order 1, but in the case of trigonometric polynomials 
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Figure 1. An example degree 30 trigonometric polynomial, along with its 1st, 
3rd, and 10th derivative. 



this is easy to see. The pth derivative of F(x) is 



N 



F^ p \x) = a n n p cos(nrr) + b n n p sin(nx), 



n=l 



where for simplicity we have assumed p is a multiple of 4, For large p the terms <2nN p cos(Nx)+ 
b N N p sin(Nx) dominate. So the zeros of F^ p \x) are close to the zeros of a N cos(Nx) + 
bjq sm.(Nx) = cjy cos(Nx + (Pn) for some real cjv and 4>Ni an d those zeros are real and equally 
spaced. 

Another general property of repeated differentiation of real entire functions of order 1 is 
that the discrepancy from equal spacing of zeros of the pth derivative scales as 0(l/p). See 
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Theorem 2.4.2 of 6J. In the case of random trigonometric polynomials we are able to obtain 
more precise information. We describe this in the next section. 

2. Statement of results 

In this section F(x) is a random trigonometric polynomial of the form (jl.l)) for which the 
ctk and bk are independent Gaussian distributed random variables with mean and identical 
variance. We wish to measure the rate at which the real zeros of the pth derivative (x) 
approach equal spacing as p — > oo. We consider the pair correlation function of the zeros of 
defined as 

(2.1) R 2;P (r) = (p p (x)p p (x + r)), 
where 

(2.2) p p (x) = 6 ( x ~ = S(F^(x))\F^ +1 \x)\. 

x k :F(p)(x k )=0 

Here 5 is the Dirac 5-function at 0, and (•) stands for expected value. Thus, p p is the density 
function of real zeros of F^ p \ and R 2jP is the density function of differences of real zeros. 

Since R 2p measures the differences between zeros, if the zeros are almost regularly spaced 
then R 2<p (x) will be large when x is close to a multiple of the average zero spacing and it will 
be small otherwise. In other words, we expect that R 2iP should approach a sum of 5-functions 
at the integers as p — ► oo. Bogomolny, Bohigas, and Lebceuf [2] obtain a general expression 
for the pair correlation function of the real zeros of a random trigonometric polynomial. 
From their results (which we describe in Section 0]) we obtain a formula for R 2>p , which we 
plot for p — 0,1, 3, 10 in Figure 121 Note that the p = case is from [2]. 

The plots in Figure 121 use the following normalization. We rescale the polynomial so 
that the average spacing between zeros is 1, that is, we are actually considering the func- 
tion F^iyix/N). As we describe in Sectional as iV — > oo the function has the expected 
fraction 

(2.3) v p -- 

of real zeros. So l/v p is the average gap between consecutive real zeros, and that is the 
spacing between the peaks in the pair correlation functions in Figure El Note that v p is also 
the density of real zeros. So the pair correlation function R 2>p (x) will equal v p on average, 
which can also be seen in Figure El 

We give an asymptotic formula for the pair correlation function R 2iP (x) as p — > oo. As 
in the general case [S| we find a 0(1 /p) discrepancy from equal spacing, and furthermore we 
find that the nearest-neighbor spacing has (appropriately rescaled) distribution function 

(2.4) 1 
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centered at 1 + m particular, the discrepancy from equal spacing is not Gaussian. The 
precise statement is 

Theorem 2.1. As p — > oo ; the pair correlation function R 2tP (x) of the real zeros of F^ 
approaches a sum of Dirac 5-functions at the nonzero integers. The 5-function near the 
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Figure 2. Plots of i? 2 , P (^). the pair correlation funtion of the zeros of the pth 
derivative F^(x), for p = 0, 1 (top row) and p = 3, 10 (bottom row). 



positive integer n is given by 
(2.5) 11 



L 2,p 



, , 1 U 
71 | 1 + — + - 

2p p 



P 



1 



+ 0(1), 



as p 



oo. 



The proof is given in Section WA\ Note that the total area under 1/(1 + 4w 2 )2 is 1, which 
shows that the above Theorem does in fact identify the ^-function near the integer n. 

Since the zeros of FW are close to equally spaced, the peak of B,2 } p(x) near x — 1 is 
almost completely due to nearest-neighbor spacings. Thus, we can read the nearest-neighbor 
distribution from the pair-correlation function, as given in (|2.4|) . In fact, we can also read 
off the next-nearest neighbor spacing (and all of the other neighbor spacings), as p — > oo, 
from the pair correlation function. Up to rescaling, all of those distributions are the same. 
This shows that there are long-term correlations between the zeros, otherwise, for example, 
the next-nearest neighbor distribution would be the convolution of the nearest neighbor 
distribution with itself. 

One motivation for this work is to understand, in general, the effect of differentiation 
on the statistics of zeros for functions which have all their zeros on a line. In particular, 
we would like to understand the effect of differentiation on the repulsion between zeros. For 
example, if f(x) has only real zeros and the zeros have the same statistics as the eigenvalues 
of the classical random matrix /3-ensemble, then the zeros of f(x) have repulsion of order j3. 
For the derivative f'{x), it is reasonable to conjecture that the zeros would have repulsion of 
order 3/5+1, because pairs of close zeros of f'(x) occur when f(x) has three closely spaced 
zeros. This topic is of interest to number theory Unfortunately, the calculations in 
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this paper do not shed light on this phenomenon because the zeros of random trigonometric 
polynomials are not in general all on the real line. We find that all derivatives of a random 
trigonometric polynomial have linear repulsion between zeros, and show that this is due to 
the fact that each derivative moves new zeros onto the real line, and those new zeros show 
linear repulsion from each other. See Section 14.21 

In the next section we discuss generalities about random trigonometric polynomials. In 
Section 0] we discuss the asymptotics of the pair correlation of the zeros of F ( p\ and in 
Section l4~2l we compute the repulsion between the zeros. 

3. Random trigonometric polynomials 

We assume that F(x) is a random trigonometric polynomial of the form ([Lip , where 
the afc and bk are independent real normally distributed random variables with mean and 
identical variance a 2 . 

The expected fraction, v p , of real zeros of F^\x) is given in (|2.3|) . This result is due to 
Dunnage [3]. It follows directly from the Kac-Rice formula, which we give in Lemma 13. 11 

Note that the original polynomial F(x) has on average 1/V3 « 57.7% real zeros, the 
4th derivative has more than 90% of its zeros real, and one must take 49 derivatives in order 
to expect 99% real zeros. Note that these are asymptotic results, and one can obtain exact 
formulas for any N. For example, the iV = 30, p — 10 example at the bottom of Figure ^ 
expects to have 96.96% real zeros, and this is approximately 1.4% larger than the asymptotic 
estimate ()2.3|) . 

We now derive ()2.3|) . 

Lemma 3.1. (The Kac-Rice Formula pHIE]j Suppose F(x) is a random trigonometric polyno- 
mial with independent real normally distributed coefficients having mean but not necessarily 
equal variance, and let 

A 2 = Var{F{x)) = (F{x) 2 ) 
B 2 = Var{F'(x)) = (F'(x) 2 ) 
C = Cov{F{x)F'{x)) = (F(x)F'(x)) 
A 2 = A 2 B 2 - C 2 . 

Then the expected number of real roots of F(x) in the interval (a, b) is 
In the case of F[x) of the form (jl.lj) we have 

N 

n=0 

JV 

B 2 =J2n 2 a 2 

n=0 

C = 0. 
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If all the coefficients of F(x) have equal variance a 2 , then the pth derivative F^ p \x) can be 
viewed as a random trigonometric polynomial of the form where the coefficients a n and 
b n have variance n 2p a 2 . So for the pth derivative we have 

N /V2P+ 1 

(3.1) A = a 2 Y,n 2p ~a 2 - 



2n + l 

n=0 y 



and 

(3.2) B = a 2 J2 



l\r2p+3 

n W „ a 2^ 

2p + 3' 

n=0 



as N — > oo. From this, formula ()2.3|) for the fraction of real roots v p follows immediately. 

4. Pair correlation of the real roots 

We use a result of Bogomolny, Bohigas, and Lebceuf |2J to compute the pair correlation 
function R 2jP of the real zeros of F^ 3 ' (x) . 

Lemma 4.1 (Appendix B of [2]). Suppose 

N 

F{x) = a n cos(nx) + b n sin(nx) 

n=0 

is a random trigonometric polynomial, where the a n and b n are independent real Gaussian 
distributed random variables with mean and variance a n . The expected value of the pair 
correlation function Rvfr) of the real zeros of F(x) is given by 



(4.1) R 2 (r) = — ^—3- (b arcsin(£ /A) + VA 2 - B 2 ) 

7T 2 C2 V ' 



where 



A = g 2 C - g x g\ 

(4.2) B = g 5 C - g 3 g 2 4 



C = g\- gl 



with 



N 



9i = ^2°l 

n=l 

N 

El 2 
n a n 

n=l 
N 

(4.3) g 3 = ^2a 2 n cos(nT) 



n=l 
N 

g A = ^no-^sin(nr) 

n=l 

95 = ^n 2 <r^cos(nr). 

n=l 



We apply the Lemma to F^^x) with a n = n p . So we have 

( 4 - 4 ) 9i ~ ^— -T and ^2 ~ TT-T^' 

2p+ 1 2p + 3 

as iV — > oo. To determine asymptotics for g 3 , g 4 and #5, we use the fact that for continuous 
functions / 

N 



n 



rl 1 / 

/ f(x)dx = lim — > /(- 

1/0 n=0 

We change variables r = nx/N, so in the variable x the mean spacing between zeros of F is 
unity. We have 

g 3 ~N 2p+1 [ cos (irxt) t 2p dt 
Jo 

(4.5) # 4 ~iV 2p+2 / sin (Trst) t 2p+l dt 

Jo 

g 5 ~iV 2p+3 / cos (irxt) t 2p+2 dt. 
Jo 

We let R,2 tP (x) denote the pair correlation function of the real zeros of the pth derivative 
F( p \x), where F(x) is given by (jl.lj) with all a,j and bj independent identical Gaussian, and 
we normalize by dividing by the square of the overall zero density N 2 /ir 2 . As iV — > 00 we 
have 

(4.6) R 2 , p {x) ~ -4 (B p arcsin(VA>) + yj Aj - 5 2 ) , 

Cp 

where 

A p = g2,pCp — 9i,p9%p 

(4.7) B p = g 5jP C p - g^ v g% v 



C p = 9 2 i, P ~ 



with 



91, P 

92, p 



1 



2p+l 
1 

2p + 3 
1 



(4.8) 93, P = / cos (nxt) t 2p dt 

Jo 

9&,p = / sin (irxt) t 2p+1 dt 
Jo 



95,) 



cos (-nxt) t 2p+2 dt. 



Plots of -R 2 , P for p = 0, 1, 3, 10 are given in Figure El 
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4.1. Large p asymptotics of the pair correlation. We determine the rate at which 
R 2 , p (x), appropriately rescaled, approaches a sum of ^-functions at the integers. 

We first find asymptotic formulas for A p , B p , and C p . Using a geometric series expansion 
for gi tP and g 2 , p , and using integration by parts followed by a geometric series expansion for 
93,p, 94,p, and g^ tP , we find (with the help of a computer algebra package), that 

A p = — (-2?rV - 2 sin(27ra;)7ra; + (4ttV - l) cos(2?nr) + l) p~ 5 + 0(p~ 6 ) 

B p = — — (cos(7ra;) + (At[ 2 x 2 — l) cos(37r:r) — 8nx sin(7nr)) p~~ 5 + 0(p~ 6 ) 

I i 

C p — - sin 2 (7ra;) p~ 2 — - ((nx cos(7rx) + sin(7r:r)) sin(7ra;)) p~ 3 



+ — (n 2 x 2 + 8 sin(27ra;)7rx + 3 (ttV - l) cos(27r:r) + 3) p~ 4 + 0(p~ 5 ) 
32 



(4.9) 



We see that A p and B p are generally of size p 5 , while C^ 2 is generally of size p 3 . But 
for x G Z we see that C p ^ 2 is of size p~ 6 . Thus, as p — * oo, if i 6 Z then R p (x) — > oo, 
otherwise R p (x) — > 0. This is exactly what one should expect because R p is approaching a 
sum of 5-functions at the integers. We now express this more precisely. 

The minima of C p are not exactly at the integers, but they are shifted over to approxi- 
mately n(l + ^) for n G Z. This is also what one would expect because the mean spacing 

of the real zeros of is approximately l/v p ~ 1 + ^. We will verify this directly from the 

above formulas. Differentiating C p with respect to x we have 

(4.10) 

C' p (x) = -— (tt {(Ap - ll)7rcos(27ra;)a; - nx + (-Ap 2 + Qp + 3ir 2 x 2 - 7) sin(27ra;))) p~ A . 

We are thinking of x as fixed and p large, so the minimum of C p is close to the solution to 

(4.11) Anpx cos(27rx) — Ap 2 sin(27ra;) = 0, 
which is equivalent to 

7TX 

(4.12) tan(27ra) = — . 

p 

The solutions with x near an integer correspond to the minima of C p , so writing x = n + £ 
and tan(27ra;) ~ 27r£, we find 

(4.13) £ ~ — , therefore x ~ n ( 1 + — 
v 1 2p \ 2p 

as expected. 

To show the shape of the ^-functions, we expand near the minima of C p . For positive 
integers n we find that 

(4.14) C p (n (l + 1 + ) = ^ (1 + 4 M 2 ) nV 4 + 0(p~ 5 ). 
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Also, 

(4.15) B p (n (l + ^ + ~J J = cos(7m)^n 2 jT 5 + 0(p- 7 ). 



Thus, arcsin(i?p/y4p) = (— l) n n/2 + 0(p ), so we obtain Theorem 12.11 

4.2. Repulsion between zeros. In j2j it was found the the real zeros of F(x) have linear 
repulsion. We find that the real zeros of (x) also show linear repulsion. 
Computing the small x asymptotics of ()4.8|) we find that 



2(2p + 3) 2 (2p + 5) (2p+l)(2p + 3)V(2p+5) 3 (2p+7) 



(4-16) ~ 



as p — > oo. So for every derivative we have linear repulsion between zeros. This appears to 
contradict the expectation that differentiation increases the repulsion between zeros. How- 
ever, there is a simple explanation. Each derivative causes more zeros to fall onto the real 
line. These "new" zeros can be closely spaced, as illustrated in Figure El 




Figure 3. Functions with zeros at | ±ia and at the integers except for 0,1. On 
the left a = 0.92, and on the right a = 1.1. Dotted curve shows the derivative. 



In Figure El all of the zeros are spaced one unit apart, except for two zeros which have 
been moved to have imaginary part ±a and real part the midpoint of the resulting gap. In 
the left plot we have a = 0.92 and in the right plot a = 1.1. Moving the complex zeros to 
have imaginary part 

(4.17) ± 2 w ±1.06975 

TX l — 8 

gives a triple zero of the derivative, so values slightly less than this give closely spaced zeros. 
Since the fraction of new real zeros for the pth derivative is 



mi, / gP + 1 /2p-l 1 

(418) V *- V *W^T3-V2FPT~ 2^ 
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which is of the same magnitude as the linear repulsion, we see that the new zeros lead to 
repulsion of magnitude x/p 2 . It seems reasonable to believe that this accounts for all of the 
linear repulsion, but we have not been able to verify this. 
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